function [E,Z2]=preparationMP(Na,E0,ph,EA,lambda,nu,f0,B,td,Kmp, EAmp, N_mp,Dmp,ha)

//1) Модель сигнала центральной антенны
exec('sglph1.sci');
for k=1:length(f0)
    f=f0(k);
   if E0(k)==0
       Ez(k,:)=zeros(1,nu);
   else
       U=sglph1(nu,td,f,B(k),ph);
       Ez(k,:)=E0(k)*U;
   end
end
//n0=1*(randn(1,nu)+1i*randn(1,nu));
n0=1*(grand(1,nu,'nor',0,1)+%i*grand(1,nu,'nor',0,1));
if length(f0)==1,
    E=Ez+n0;
else
    E=sum(Ez,'r')+n0;  
end;

 // 2) Модель сигнала системы периферийных антенн
L =0.5* lambda;  
// L=0.2*lambda;  
H=0;
// en=randn(Na,nu)+%i*randn(Na,nu); 
//en=1*(randn(Na,nu)+1i*randn(Na,nu));
en=1*(grand(Na,nu,'nor',0,1)+%i*grand(Na,nu,'nor',0,1));

exec('Ant_ran.sci');
[Z1,Z2]=Ant_ran(Ez,en,Na,L,EA,lambda,H);

// [Ez, Z2]=ant_ran_mp_ground(Ez,Z2,Na,L,lambda,H,EA, td,ha) ;

if  N_mp>0
    exec('ant_ran_mp2.sci');
    [Ez, Z2]= ant_ran_mp2(Ez,Z2,Na,L,lambda,H,Kmp,Dmp,EA, EAmp, N_mp,td);
end;
tau_max=5000;
E(1:tau_max)=[];
Z2(:,1:tau_max)=[];

endfunction
